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TITLE OF THE INVENTION 

A RECEIVING APPARATUS , A RECEIVING METHOD, 
AND A RADIO COMMUNICATIONS SYSTEM 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a 
receiving apparatus, a receiving method, and a radio 
communications system for demodulating (separating) 
a signal using two or more antennas. 

2 . Description of the Related Art 
In the fourth generation mobile 

communications, a radio communications method that 
is capable of providing a high transmission speed is 
required. From this viewpoint, a MIMO (multiple 
inputs and multiple outputs, i.e., two or more I/O) 
multiplexing method attracts attention, wherein each 
of antennas (for transmission and reception) 
transmits a. different signal from others at the same 
time and the same frequency using MIMO channels. 

The MIMO multiplexing method is described 
with reference to Fig. 14 that shows the 
configuration of a MIMO communication system using 
two or more antennas. According to this system, the 
transmission speed is increased in proportion to the 
number of transmitting antennas without increasing a 
transmission bandwidth by antennas lOlli through 
1011 N transmitting respective different signals at 
the same frequency, and antennas 1021i through 1021 N 
simultaneously receiving all the different signals 
(for example, Non-patent Reference 1). 

Further, BLAST (Bell Labs Layered Space - 
Time) is one of the technologies generically named 
MIMO. According to BLAST, different signals are 
simultaneously transmitted at the same frequency, 
i.e., a parallel transmission, from two or more 
transmitting antennas; and, on the receiving side, 
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the signals are separated by diversity reception 
that is controlled by interference suppression and 
by a replicas subtraction (for example, Non-patent 
Reference 1) . 

5 As described above, while it is possible 

to realize a high transmission speed by the MIMO 
multiplexing method, since different data sequences 
are transmitted from the transmitting antennas at 
the same frequency and at the same time slot, a 

10 receiving unit is required to separate the signals 
in order to extract the data sequences transmitted 
from each transmitting antenna from the received 
signal for demodulation . 

Various methods are proposed for 

15 separating the MIMO multiplexed signals. For example, 
according to a signal separation algorithm using a 
linear filter such a minimum mean square error 
(MMSE) method, and a zero forcing (ZF) method, 
signals received by antennas, the number of which 

20 antennas is equal to or greater than the number of 
the transmitting antennas , are compounded, such that 
receiving power from transmitting antennas other 
than a target transmitting antenna is suppressed 
(minimized) . According to this method, the amount of 

25 operations on the receiving side is relatively small. 

Further, according to a signal separation 
algorithm using a maximum likelihood detection 
method (MLD) , replica candidates of signals received 
from each transmitting antenna are generated, and a 

30 replica of a received signal that gives the smallest 
Euclidean distance between the received signal and a 
sum of the replica candidates of the signals 
received from all the transmitting antennas is 
obtained. In this way, the MLD method gives the most 

35 probable signal sequence of each transmitting 
antenna. The MLD method provides more accurate 
signal separation and superior demodulation 



performance as compared with the MMSE method; 
however, the amount of operations required of signal 
separation exponentially increases as the number of 
the antennas is increased. Then, a proposal (for 
example, Non-patent Reference 2) is made wherein the 
amount of operations is decreased in the MLD by 
greatly decreasing the number of signal point 
candidates for calculating the square Euclidean 
distance using QR factorization. 

By the way, conventional direct spreading 
(DS) CDMA is a communication method wherein a 
secondary modulation is performed for spectrum - 
spreading a signal that is modulated by conventional 
information data with a high - speed spreading 
signal such that two or more parties can communicate 
using the same frequency band. In the radio 
communications using DS-CDMA, as the bandwidth used 
for signal transmission becomes great, multipath 
fading (frequency selective fading) is generated, 
and a transmission signal is received through two or 
more paths (multipath) having different propagation 
delay times . 

While receiving quality can.be improved in 
DS-CDMA with a rake reception method, wherein two or 
more signals through the multipath are compounded, 
interference (multipath interference) is generated 
between different paths. For this reason, the 
receiving quality improvement by the rake reception 
method is offset by the multipath interference. 

The magnitude of the multipath 
interference is proportional to the inverse number 
of a spreading rate that is defined by a ratio of a 
chip rate (representing a speed of multiplication by 
the spreading signal) to a symbol rate of 
information symbols. For this reason, if the 
spreading rate is brought close to 1 in an attempt 
to increase an information bit rate, degradation of 
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the receiving quality due to the multipath 
interference becomes dominant rather than the 
improvement obtained by the rake reception method. 
This poses a problem in that the receiving quality 
5 is degraded when transmitting at a high speed. Then, 
in an attempt to cope with the problem, a multipath 
interference canceller is proposed (for example, 
Non-patent Reference 3) . 

According to the multipath interference 

10 canceller as disclosed by Non-patent Reference 3, 

high quality reception under a multipath environment 
is realized by 

estimating a received signal sequence for 
every path based on an signal sequence that is 

15 estimated by. a provisional rake receiving result and 
a channel coefficient (complex envelope of a 
propagation path) of each reception path, 

deducting all estimated signal sequences 
obtained from paths other than a certain path (a 

20 target path) from the received signal, and 

repeating the deducting operation for all 
the paths. In this way, the multipath interference 
of the signal of every path is reduced. Then, using 
the signals, for the multipath interference of which 

25 signals have been reduced, a final rake reception is 
carried out . 

Further, 2-dimensional MMSE is proposed as 
a signal separating method, whereby the multipath 
interference in MIMO multiplexing in the radio 

30 communications using DS-CDMA is reduced, With the 2 — 
dimensional MMSE, both interference generated by 
received signals from other transmitting antennas 
and multipath interference are simultaneously 
reduced . 

35 [Non-patent Reference 1] G. J. Foschini, 

Layered Space-Time Architecture for Wireless 
Communication in a Fading Environment When Using 
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Multiple Antennas , Bell Labs Technical Journal, Vol. 
1, No. 2, autumn 1996, pp 41 through 59. 

[Non-patent Reference 2] Bin Dong, Xaodong 
Wang , and Amaud Douce t , "Sampling- based Near- 
5 optimal MIMO demodulation Algorithms" in Proc. 42nd 
IEEE Conference on Decision and Control, Hawaii, 
Dec . , 2003 

[Non-patent Reference 3] K. Higuchi, A. 
Fujiwara, and M . Sawahashi , "Multipath Interference 

10 Canceller for High-Speed Packet Transmission With 
Adaptive Modulation and Coding Scheme in W-CDMA 
Forward Link, M IEEE J. Select. Areas Commun., Vol.20, 
No. 2, pp.419 through 432, February 2002. 

[Non-patent Reference 4] Frederik Petre et. 

15 al , "Combined Space-Time Chip Equalization And 
Parallel Interference Cancellation For DS-CDMA 
Downlink With Spatial Multiplexing" in Proc. IEEE 
PIMRC 2002 . 

2 0 SUMMARY OF THE INVENTION 

[Problem (s) to be solved by the Invention] 
When applying the MIMO multiplexing method 
as described above for increasing the information 
bit rate of the radio communications using DS-CDMA, 
25 the spreading rate should be made small such that 
the information bit rate per transmitting antenna 
can be increased. In this case, as shown in Fig. 15, 
there are generated 

mutual interference between the 
30 transmitting antennas of the paths of the same 

receiving timing (indicated by (a) in Fig. 15) , and 

the multipath interference (indicated by 
(b) and (c) in Fig. 15) from the multiple paths of 
all the transmitting antennas of different receiving 
35 timings. For this reason, accuracy of the signal 

separation on the receiving side is greatly degraded 
due to the multipath interference as in the case of 



the rake reception with one antenna transmission 
described above. 

The proposed multipath interference 
canceller described above is configured in 
consideration of a signal transmitted from only one 
antenna. If this is applied to MIMO multiplexing, 
accurate estimation of a signal sequence with the 
rake reception cannot be obtained. 

Further, as evident from Fig. -15, the 
multipath interference in the case of MIMO 
multiplexing is generated not only by one antenna, 
but also by two or more antennas. Accordingly, it is 
necessary to reduce the multipath interference due 
to multiple antennas for highly accurate reception. 

However, the signal separation of the MMSE 
method theoretically is poorer than the MLD method, 
and a desired performance cannot be obtained on the 
receiving side. Further, according to the 2- 
dimensional MMSE method, a filter coefficient is 
controlled in an attempt to reduce both the 
interference generated by the received signals from 
other transmitting antennas and the multipath 
interference. For this reason, signal separation 
accuracy is further degraded, posing a problem. 

The present invention is made in view of 
the problems described above, and offers a receiving 
apparatus, a receiving method, and a radio 
communications system that are capable of realizing 
highly accurate signal separation by reducing 
degradation of the receiving quality due to the 
multipath interference . 

[Means for solving the Problem} 

According to an embodiment of the present 
invention, a receiving apparatus is provided for 
receiving a signal using the CDMA method, wherein 
signals are transmitted from M transmitting antennas 
(M is a positive integer) , and received by N 



receiving antennas (N is a positive integer) . The 
receiving apparatus includes: 

multipath receiving signal demodulating 
units configured to 

demodulate (primary demodulation) the 
signals received by the receiving antennas, 

estimate the transmitted signals from 
the transmitting antennas , and 

obtain a received signal of each path 
of the receiving antennas in a multipath 
environment based on the estimation; 

multipath interference canceling units 
configured to deduct the obtained received signals 
of paths other than a target path from the signals 
received by the receiving antennas to obtain a 
multipath interference cancelled signal; and 

a demodulator configured to perform a 
secondary demodulation of the multipath interference 
cancelled signal output by the multipath 
interference canceling unit. 

According to an aspect of the present 
invention, the multipath receiving signal 
demodulating unit of the receiving apparatus carries 
out the primary demodulation using a minimum mean 
square error (MMSE) method. 

According to an aspect of the present 
invention, the multipath receiving signal 
demodulating unit of the receiving apparatus carries 
out the primary demodulation using a maximum 
likelihood detection (MLD) method. 

According to an aspect of the present 
invention, the multipath receiving signal 
demodulating unit of the receiving apparatus carries 
out the primary demodulation of two or more paths in 
one block using the maximum likelihood detection 
method that uses QR factorization. 

According to an aspect of the present 



invention, the multipath receiving signal 
demodulating unit of the receiving apparatus carries 
out the primary demodulation of each of the paths 
using the maximum likelihood detection method that 
uses QR factorization . 

According to an aspect of the present 
invention, the multipath receiving signal 
demodulating unit of the receiving apparatus 
controls the amplitude of the received signal based 
on the correctness 1 ikelihood of a transmission 
symbol sequence estimated using any one of the 
methods described above. 

According to an aspect of the present 
invention, the multipath receiving signal 
demodulating unit of the receiving apparatus 
estimates a channel coefficient using a known pilot 
signal transmitted from the M transmitting antennas. 

According to an aspect of the present 
invention, the receiving apparatus includes a 
predetermined number of multipath receiving s ignal 
demodulating units, and a predetermined number of 
multipath interference canceling units connected in 
stages (cascaded) . 

According to an aspect of the present 
invention, when the receiving apparatus includes the 
staged connection of the multipath receiving signal 
demodulating units, the estimated channel 
coefficient is updated in each stage based on the 
known pilot signal transmitted from the M 
transmitting antennas using the multipath 
interference cancelled signal. 

Further, according to an aspect of the 
present invention, the demodulating unit of the 
receiving apparatus carries out the secondary 
demodulation using the maximum likelihood detection 
method . 

According to an aspect of the present 
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invention, the demodulating unit of the receiving 
apparatus carries out the secondary demodulation of 
two or more paths in one block using the maximum 
likelihood detection method that uses QR 
5 factorization. 

According to an aspect of the present 
invention, the demodulating unit of the receiving 
apparatus carries out the secondary demodulation of 
each of the paths using the maximum likelihood 
10 detection method that uses QR factorization. 

According to an aspect of the present 

invention, 

when code multiplexed signals are 
transmitted from the M transmitting antennas, the 

15 multipath receiving signal demodulating unit of the 
receiving apparatus carries out the primary 
demodulation of a signal received by each receiving 
antenna, and the received signal of each path for 
every receiving antenna is obtained for every 

20 spreading signal, 

the multipath interference canceling unit 
generates a signal by deducting the obtained signal 
corresponding to all the spreading signals of paths 
other than a target path from the signals received 

25 by each of the receiving antennas to obtain a 
multipath interference cancelled signal, and 

the demodulating unit carries out a 
secondary demodulation of the multipath interference 
cancelled signal for every spreading signal. 

30 [Effect of the Invention] 

According to the embodiment of the present 
invention, when different data are simultaneously 
transmitted from two or more transmitting antennas 
using the CDMA method, highly accurate separation of 

35 the signals transmitted from different transmitting 
antennas, the multipath interference of which 
signals is reduced, can be realized. Consequently, 
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the receiving quality in the multipath fading 
environment is greatly improved. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 Fig. 1 is a block diagram of a radio 

communications system including a receiving 
apparatus according to an embodiment of the present 
invention. 

Fig. 2 is a block diagram of a first 
10 embodiment of a multipath receiving signal 
demodulator . 

Fig. 3 is a block diagram of a second 
embodiment of the multipath receiving signal 
demodulator . 

15 Fig. 4 is a block diagram of a third 

embodiment of the multipath receiving signal 
demodulator . 

Fig. 5 is a block diagram of a fourth 
embodiment of the multipath receiving signal 
20 demodulator. 

Fig. 6 is a block diagram of a channel 
coefficient estimating unit. 

Fig. 7A is a data map showing an example 
of a frame structure of a signal transmitted from a 
25 transmitting apparatus. 

Fig. 7B is a data map showing another 
example of the frame structure of the signal 
transmitted from the transmitting apparatus. 

Fig. 8 is a block diagram of a second 
30 embodiment of the receiving apparatus according to 
the present invention. 

Fig. 9 is a block diagram of a first 
embodiment of a demodulator according to the present 
invention . 

35 Fig. 10 is a block diagram of a second 

embodiment of the demodulator according to the 
present invention . 
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Fig. 11 is a block diagram of a third 
embodiment of the demodulator according to the 
present invention. 

Fig. 12 is a block diagram of a third 
5 embodiment of the receiving apparatus according to 
the present invention. 

Fig. 13 is a graph showing results of a 
computer simulation about the present invention. 

Fig. 14 is a block diagram of a MIMO 
10 communication system using two or more antennas. 

Fig. 15 is a schematic diagram showing 
that accuracy of demodulation (signal separation) 
deteriorates due to multipath interference in the 
case of MIMO multiplexing of DS-CDMA. 

15 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

. [Description of Notations] 
1 Radio Communications System 
10 1010 Transmitting apparatus 
20 Hi through ll m , lOlli through 1011 n 

Transmitting antenna 

20, 200, 600, 1020 Receiving apparatus 
21i through 21 n , 1021i through 1021 n 
Receiving antenna 
25 22, 30 through 50, 70, 211 through 213, 

611 through 614 Multipath receiving signal 
demodulator 

23, 221i, 221 2 , 222i, 222 2 , 223i, 223 2 , 615 
through 618 Multipath interference canceling unit 
30 24, 300, 400, 500, 621, 622 Demodulator 

31, 41, 51, 71, 100, 311, 411, 511 Channel 
coefficient estimating unit 

32 Linear filter coefficient Calculating 

unit 

35 33 Linear filter 

34i through 34 m , 42 through 45, 52 through 
55, 72 through 75, 312 through 315, 412 through 415, 



- 12 - 



512 through 515 Reverse spreading unit 

35i through 35 m Transmission symbol 

sequence estimating unit 

36i through 36 m , 50i through 50 m , 62i 
5 through 62 m , 84i through 84 m Multipath receiving 

signal demodulator 

46, 58, 80, 316, 418, 520 Transmission 
symbol candidate generating unit 

47, 317 Received signal replica generating 

10 unit 

48, 60, 82, 318, 420, 522 Likelihood 
calculating unit 

49i through 4 9 m , 61i through 61 ra , 83i 
through 83 m Transmission symbol sequence estimating 
15 unit 

56, 76, 77, 416, 516, 517 QR factorizing 

unit 

57, 78, 79, 417, 518, 519 Q H operating 

unit 

20 59, 81, 419, 521 Conversion signal replica 

generating unit 

101 through 104 Correlation detecting unit 
111 through 114 Replicas signal generating 

unit 

25 319, 421, 523 Transmitting sequence 

.estimating unit 

631 through 634, 641 through 644 Adder 
[Best Mode of Carrying Out the Invention] 
Hereafter, embodiments of the present 
30 invention are described with reference to the 
attached drawings . 

Fig. 1 is a block diagram of a radio 
communications system 1 including a receiving 
apparatus 20 according to the embodiment of the 
35 present invention. The radio communications system 1 
according to the present invention uses one of a 
MIMO channel and a MISO (Multi Input Single Output, 



i.e., two or more inputs and a single output, using 
only one receiving antenna) channel in DS-CDMA . In 
the following, an example wherein a signal is 
transmitted using a MIMO channel is described. 

As shown in Fig. 1, the radio 
communications system 1 includes a transmitting 
apparatus 10 and the receiving apparatus 20 that are 
connected through radio communications. The 
transmitting apparatus 10 carries out serial to 
parallel conversion of a data bit sequence that is 
input for transmission into M transmitting sequences 
modulates the transmitting sequences (data 
modulation) using the same frequency band and the 
same spreading signal, and simultaneously transmits 
the modulated sequences from M transmitting antennas 
Hi through ll m as transmission signals. Each of the 
transmission signals is propagated through L paths 
(multipath) , and is received by N receiving antennas 
21i through 21 n that are connected to the receiving 
apparatus 20. The receiving apparatus 20 includes a 
multipath receiving signal demodulator 22, N 
multipath interference canceling units 23i through 
23 n , and a demodulator 24. 

Here , a trans mis sion signal from a 
transmitting antenna ll m connected to the 
transmitting apparatus 10 is expressed by S m (t) , and 
a received signal received by a receiving antenna 
21 n connected to the receiving apparatus 20 is 
expressed by r n (t). Then, r n (t) can be described as 
follows . 

[Expression 1] 

Here, h m , n ,<(t) is a channel coefficient of 
a reception path I when receiving a transmission 
signal from the transmitting antenna ll m by the 



receiving antenna 21 n , and t< is a propagation delay 
time of the reception path I . 

Next, operations of the receiving 
apparatus 20 according to the present invention are 
described . 

r 

With the receiving apparatus 20, signals 
r n (t) received by the receiving antennas 21i through 
21 n are provided to the multipath receiving signal 
demodulator 22. The multipath receiving signal 
demodulator 22 carries out provisional signal 
separation (a primary demodulation) based on the 
signals received by the N receiving antennas 21i 
through 21 n such that received signal sequences as 
described by the following Expression 2 for the 
reception paths between the transmitting antennas 
Hi through ll m and the receiving antennas 21i 
through 21 n are estimated and output. 

[Expression 2 ] 

An algorithm for carrying out the primary 
demodulation is described below. 

Next, the signals received by the 
receiving antennas 21i through 21 n , and the received 
signal sequences of every reception path from the 
transmitting antennas Hi through ll m are provided to 
the corresponding multipath interference canceling 
units 23i through 23 n . The multipath interference 
canceling units 23i through 23 n (that number N, the 
same number as the receiving antennas 21i through 
21 n ) obtain received signals after the multipath 
interference cancellation r n ,e(t) by deducting all 
the received signals of other paths than a target 
path from the received signal according to the 
following expression, and output the multipath 
interference cancelled signals r n ,<(t) . 
[Expression 3] 



The demodulator 24 receives NxL multipath 
interference cancelled signals r n ,i(t) provided by 
the multipath interference canceling units 23i 
through 23 n . The demodulator 24 outputs a secondary 
demodulation result that is a likelihood , the 
likelihood being one of a likelihood of the 
transmission symbol sequence and a likelihood of a 
bit transmitted by the transmission symbol sequence. 
The secondary demodulation result is provided to an 
error correction (channel) decoder, such as a 
Viterbi decoder and a turbo decoder, for decoding. 

As described above, according to the 
present embodiment, the multipath receiving signal 
demodulator 22 carries out the primary demodulation, 
and the transmission signals are estimated. Then, 
the received signal for every path is estimated by 
multiplying the estimated received signal and the 
channel fluctuation values (channel coefficients) . 
Then, the multipath interference canceling units 23 
subtract estimated received signals other than the 
signal received through the target path. In this way, 
the demodulator 24 is able to demodulate the 
received signal after removing multipath 
interference, and signal separation can be performed 
with high accuracy. 

That is, even if the DS— CDMA method is 
applied to an uplink, and the MIMO multiplexing 
method is used, degradation of the signal separation 
accuracy resulting from the multipath interference 
can be avoided. 

Fig. 2 is a block diagram showing a 
multipath receiving signal demodulator 30 according 
to a first embodiment of the multipath receiving 
signal demodulator 20 shown in Fig. 1. According to 
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the first embodiment, the multipath receiving signal 
demodulator 30 uses the MMSE algorithm as the 
primary demodulation method. 

As shown in Fig. 2, the multipath 
5 receiving signal demodulator 30 includes a channel 
coefficient estimating unit 31, a linear filter 
coefficient calculating unit 32, a linear filter 33, 
M reverse spreading units 34i through 34 m , M 
transmission symbol sequence estimating units 35i 

10 through 35 m , and M multipath receiving signal 
demodulators 36i through 36 m . 

According to the multipath receiving 
signal demodulator 30 of the first embodiment, 
signals r n (t) received by N receiving antennas 21i 

15 through 21 n are provided to the channel coefficient 
estimating unit 31, and channel coefficients h m , n ,c 
of paths between the receiving antennas 21i through 
21 n and the transmitting antennas Hi through ll n 
(Fig. 1) are estimated. 

20 Then,, the linear filter coefficient 

calculating unit 32 calculates a coefficient for the 
linear filter 33, which carries out a filtering 
process for equalization using the linear filter 
-coefficient. The process described above can be 

25 realized by signal processing in one of the time 

domain and the frequency domain; in the following, 
descriptions are presented about an example wherein 
signal processing is carried out in the frequency 
domain . 

30 The linear filter coefficient calculating 

unit 32 calculates the linear filter coefficient 
that simultaneously minimizes delay path components 
of signals from a desired transmitting antenna and 
receiving signal components from other transmitting 

A 

35 antennas using the estimated channel coefficient h m ^ 
as described by Expression 4. 



[Expression 4] 

A 

H m ,n J 

Specifically, an impulse response of a 
channel (channel impulse response) between the 
transmitting antennas Hi through ll m and the 
receiving antennas 21i through 21 n is calculated 
using the estimated channel coefficient as described 

A 

by Expression 5 h m *j estimated using a pilot channel 
in the time domain, and a delay time of each 

A 

reception path r, as described by Expression 6. 
[Expression 5] 

A 

hmj*J 

[Expression 6] 



Next, FFT for N f points is carried out on 
the channel impulse response obtained as described 
above. Here, N f is a number that is equal to a 
product of 

the number of chips equivalent to a block 
size that are linear-filtered, and the number of 
times of over-sampling. In this way, the frequency 

components as described by Expression 7 of the 

channel coefficients of the transmission signals 
from the transmitting antennas Hi through ll m to 
the receiving antennas 21i through 21 n are estimated 
[Expression 7] 

run 

Using this, an estimated value of a' 
channel matrix having N rows x M columns as 
described by Expression 9 is obtained as described 
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by Expression 8. 



[Expression 8] 



H </) =[h^](l<n<N,l<m<M) 



[Expression 9] 



5 



Using Expression 10, linear filtering 
coefficients are calculated for every frequency 
component after FFT . 

[Expression 10] 



10 



The linear filter coefficients can be 
determined in accordance with one of ZF criterion 
and MMSE criterion. 

For example, the linear filter 
15 coefficients by the ZF criterion are obtained 
according to the following expression. 



Further, the linear filter coefficients by 
20 the MMSE criterion are obtained according to the 
following expression . 

[Expression 12] 



Here, the last term of Expression 12 that 
25 is described by Expression 13 represents a noise 
component . 

[Expression 13] 



30 of N sequence received signals into a received 



[Expression 11] 





Next, the linear filter 33 converts each 



signal Y 



(f ) 



in the frequency domain by carrying out 
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10 



15 



20 



25 



FFT for the N f points. Then, Y (f> is multiplied by 
W (f> for obtaining estimated values as described by 
Expression 14 of M transmission signals according to 
Expression 15, where both 

collapse of the orthogonality between code 
channels due to channel fluctuation in the frequency 
domain (MPII (Multi-Path Interference) in the time 
domain) , and 

interference between the transmitting 



antennas 



are simultaneously equalized (reduced) . 
[Expression 14] 

[Expression 15] 
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In addition 
[Expression 16] 

is reconverted into every provisional 
demodulating signal (primary demodulation signal) S m 

as described by Expression 17 in the time domain by 
carrying out parallel-to-serial conversion after 
IFFT for the N f points. 

[Expression 17] 



The acquired provisional demodulating 
sequences (primary demodulation sequences) for the 
transmitting antennas are provided to the 
corresponding reverse spreading units 34i through 
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34 m . The reverse spreading units 34i through 34 m 
carry out reverse spreading of the provisional 
demodulating sequences for the transmitting antennas 
with the same spreading signal as used at the time 
5 of transmission, and reverse-spread signals z m are 
output to the corresponding transmission symbol 
sequence estimating units 35i through 35 m . The 
transmission symbol sequence estimating units 35i 
through 35 m carry out one of a hard decision and a* 
10 soft decision of the corresponding reverse-spread 
signals z m such that transmission symbol sequences 
as described by Expression 18 are estimated and 
output. 



The transmission symbol sequence 
estimating units 35i through 35 m may carry out the 
hard decision as follows. 

Here, a transmission symbol is expressed 
20 by d m ,i, where "i" is a symbol candidate number 

satisfying Expression 19 (1</<C), and C expresses 
the number of transmission symbol points, that is, 4 
for QPSK and 16 for 1 6QAM . 

[Expression 19] 



A symbol candidate number "i" that makes 
the transmission symbol d m ,i closest to z m (refer to 
Expression 20) is determined. 

[Expression 20] 



[Expression 18] 



15 



25 



(1</<C) 



30 




Using the "i" determined as above, a 
determination of Expression 21 is made. 
[Expression 21] 




= d 
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Alternatively, the transmission symbol 
sequence estimating units 35i through 35 m may carry 
out the soft decision as follows. 

First, a bit sequence for the soft 
decision is obtained according to Expression 22. 

[Expression 22] 



u mj = tanh 



mJ 



Here, the expression A mJ shown by 

Expression 23 is a logarithmic likelihood ratio of a 
10 bit "i" of the signal from the transmitting antenna 

Urn. 

[Expression 23] 



15 



20 



The logarithmic likelihood ratio 



A mJ can 



be described by Expression 24 
[Expression 24] 



A =ln 



P(b mJ = +l\z m ) 



P(b mJ =-l\z m ) 
1 



2cr 2 (n) 



Here, S min ,v expresses a symbol candidate 
having the smallest Euclidean distance from the 
transmission signal point z m out of a set of symbols 
whose i-th bit is '• v" , and a 2 is noise power. 

Then, a soft decision symbol 

[Expression 25]. 



25 is estimated as shown by Expression 27 

using the expression of Expression 26. 
[Expression 26] 
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u 



20 



[Expression 27] 
Providing : d m = x m +j -y m 



IncaseofQPSK < 



In case of 1 6QAM 



Finally, the multipath receiving signal 
5 demodulators 36i through 36 m receive the 

transmission symbol sequences as described by 

Expression 28 d m that are estimated as described 

above. 

[Expression 28] 

10 d m 

Then, the input transmission symbol 

sequences as described by Expression 29 d m are 

multiplied by the spreading signal and the channel 
coefficients such that a received signal sequence 

15 {as shown by Expression 30} { I m „jQ) ) for every 

reception path from the transmitting antennas to the 
receiving antennas is estimated as shown by 
Expression 31, and output. 

[Expression 29] 



[Expression 30] 
[Expression 31] 
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Next, another embodiment of the multipath 
receiving signal demodulator is described. 

Fig. 3 is a block diagram showing a 
5 multipath receiving signal demodulator 40 that is a 
second embodiment of the multipath receiving signal 
demodulator of the present invention. According to 
the second embodiment, an MLD algorithm is used for 
the multipath receiving signal demodulator 40 as the 
10 method of the primary demodulation. 

The multipath receiving signal demodulator 
40 includes a channel coefficient estimating unit 41, 
NxL reverse spreading units 42 through 45, a 
transmission symbol candidate generating unit 46, a 
15 received signal replica generating unit 47, a 
likelihood calculating unit 48, M transmission 
symbol sequence estimating units 49i through 49 m , 
and M multipath receiving signal demodulators 50i 
through 50 m . 

20 According to the multipath receiving 

signal demodulator 40 of the second embodiment , 
signals r n (t) received by the receiving antennas 21i 
through 21 n are provided to the channel coefficient 
estimating unit 41, and the channel coefficients 

25 h m , n ,* of the paths between the receiving antennas 21i 
-through 21 n and the transmitting antennas Hi 
through ll n are estimated. 

Next, the reverse spreading units 42 
through 45, which number NxL, perform reverse 

30 spreading of the received signals r n (t) with the 

same spreading signal as the spreading signal used 
at the time of transmission at the receiving timing 
of each path, and NxL reverse-spread signals z n ,f for 
the paths of the receiving antennas are obtained. 

35 The transmission symbol candidate 

generating unit 46 generates and outputs 



transmission symbols d m#i corresponding to each 
transmitting antenna, where "i" is a symbol 
candidate number that suffices for Expression 32. 

[Expression 3 2 ] 

(1</<C) 

where C is the number of transmission 
symbol points, for example, C=4 for QPSK, and C=16 
for 16QAM. 

The received signal replica generating 
unit 47 receives the transmission symbols d m ,i 
generated by the transmission symbol candidate 
generating unit 46, and the channel coefficients 
estimated by the channel coefficient estimating unit 
41 such that received signal replicas as described 
by Expression 33 are generated according to 
Expression 34 and output. 

[Expression 33] 

z 

[Expression 34] 

A A A 

= K: d XA +K„; +^,/ d M4i 

The likelihood calculating unit 48 
receives the reverse-spread signals z n ,t and the 
received signal replicas as described by Expression 
35 . 

[Expression 35] 

Z , . . . 

Then, the likelihood calculating unit 48 
calculates errors according to Expression 36. 
[Expression 36] 



The transmission symbol sequence 
estimating units 49i through 49 m receive the 
generated transmission symbols corresponding to the 



transmitting antennas, receive the corresponding 
error signals, determine the minimum of the errors, 
and estimate transmission symbol sequences, as 
described by Expression 37, which give the minimum 
errors . 

[Expression 37] 



In the case where the transmission symbol 
sequence estimating units 49i through 49 m carry out 
hard decisions, Expression 38 is used. 

[Expression 38] 

d m = d m ,un where h > h > * '« give minimum 

When the transmission symbol sequence 
estimating units 49i through 49 m carry out soft 
decisions, the following process can be used. 

First, a bit sequence for soft decision is 
obtained according to the following expression . 

[Expression 39] 

f 2 \ 



u 



= tanh 



. 2 



Here , 

[Expression 40] 



is a logarithmic likelihood ratio of a bit 
"i" of a transmitting antenna ll m , and is expressed 
as follows . 

[Expression 41] 



" P(b mJ =-l\z ul ,.., Ztl J 

~ I ( e -a ) 

2<T 7 (n) 

Here, e m i n , v is the minimum value of 
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[Expression 42] 



n.*2 

wherein the i-th bit is 1 
expresses noise power. 

The soft decision symbol 
[Expression 43] 



and cr (n) 



10 



15 



20 



is estimated as shown by Expression 45 
using the expression of Expression 44. 
[Expression 44] 

[Expression 45] 
Providing: d m = £ m +y -y m 

In the case of QPSK J ^ m ~~ ^ 

[y 



= u 



m,2 



In the case of 16QAM 



) 



Finally, the multipath receiving signal 
demodulators 36i through 36 m (50i through 50 m ) 
receive the transmission symbol sequences described 
by Expression 46 that are estimated as described 
above . 

[Expression 46] 



Then, the received transmission symbol 
sequences, Expression 47, are multiplied by the 
spreading signal and the channel coefficients. 

[Expression 47] 



25 



In this way, the received signal sequences 
described by Expression 48 for every reception path 
of the receiving antennas from the transmitting 
antennas are estimated according to Expression 49, 
and output. 

[Expression 48] 

Lao 

[Expression 49] 

. Fig. 4 is a block diagram of a multipath 
receiving signal demodulator 50 that is a third 
embodiment of the multipath receiving signal 
demodulator. According to the third embodiment, the 
MLD algorithm is used, where the multipath receiving 
signal demodulator 50 carries out the primary 
demodulation of the paths in one block. 

As shown in Fig. 4, the multipath 
receiving signal demodulator 50 includes a channel 
coefficient estimating unit 51, NxL reverse 
spreading units 52 through 55, a QR factorizing unit 
56, a Q H operating unit 57, a transmission symbol 
candidate generating unit 58, a conversion signal 
replica generating unit 59, a likelihood calculating 
unit 60, M transmission symbol sequence estimating 
units 61i through 61 m , and M multipath receiving 
signal demodulators 62 1 through 62 m . 

According to the multipath receiving 
signal demodulator 50 of the third embodiment, the 
channel coefficient estimating unit 51 receives the 
signals r n (t) received by the receiving antennas 21i 
through 21 n , and the channel coefficients h m , n ,< of 
the paths between the receiving antennas 21i through 
21 n and the transmitting antennas Hi through ll n are 
estimated . 

Next, the reverse spreading units 52 
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15 



20 



through 55 , which number NxL # perform reverse 
spreading of the received signals r n (t) with the 
same spreading signal as used at the time of 
transmission at the receiving timing of each path, 
and NxL reverse-spread signals z n ,G of the paths of 
the receiving antennas are obtained. 

Next, the QR factorizing unit 56 generates 
a channel matrix as shown by Expression 50, the 
channel matrix including channel coefficients 
expressed in NxL rows x M columns, performs QR 
factorization of the channel matrix, and a Q matrix 
and an R matrix are output. 

[Expression 50] 



Channel Matrix H = 



^1.1,1 


*2.,.l 


*3.1.1 


^4.1.1 


^1.2.1 


*2.2,, 


^3.2.1 


*4.2.1 


^1,3.1 


^2.3.1 


^3.3.1 


*4,3.1 


^1.4.1 


^2.4,1 


^3.4.1 


^4,4.1 


^1.1.2 


^2.1,2 


^3,1.2 


^4,1.2 


^1.2.2 


^2.2.2 


^3.2.2 


^4,2,2 


^1.3,2 


^2.3.2 


^3.3.2 


^4,3.2 


^1.4.2 


^2,4,2 


^3.4.2 


^4.4.2 


4, N 


= 4, L 


= 2) 





( Example where N 
H = QR 

The Q matrix output from the QR 
factorizing unit 56 is a unitary matrix having (NxL) 
rows x M columns, and suffices for Q H Q=I . 

Here, H expresses conjugate complex 
transposition, and "I" expresses a unit matrix. 
Further, the R matrix turns into an upper triangular 
matrix having M rows x M columns. 

The Q H operating unit 57 performs 
operations according to the following expression. 

[Expression 51] 
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*1 




X = 




= Q H Z = 










*4_ 





9u 



01.2.1 01.3.1 01.4.1 01.1,2 01.2.2 01.3,2 01.4.2 



02.1.1 02.2.1 02.3.1 02.4.1 02.1.2 02.2.2 02.3.2 02.4.2 

• • * • 
03.1,2 03,2,2 03.3.2 03.4,2 

• • * • 
04.1,2 04.2.2 04.3.2 04.4.2 



03.L1 03.2.1 03.3.1 03.4.1 

* * * • 

04.1.1 04.2,1 04.3.1 04.4.1 



Q H (HZ) + N) = Q H (QRD + AD = RD + Q H N = 



Example where M=4, N=4, L = 2 



*4.1 
*1.2 
*2.2 
*3.2 
*4.2 





r i2 


>13 


*14 


0 




^23 


>24 


0 


0 




>34 


0 


0 


0 


r*4 



d 2 



10 



15 



The transmission symbol candidate 
generating unit 58 generates the transmission 
symbols d M/ i for the M transmitting antennas 11, and 
outputs them to the conversion signal replica 
generating unit 59. The conversion signal replica 
generating unit 59 receives 

the transmission symbols generated by the 
transmission symbol candidate generating unit 58 , 
and 

the R matrix provided by the QR 
factorizing unit 56 

for generating conversion signal replicas 
as described by Expression 52. 

[Expression 52] 



M J 

The conversion signal replicas are 
generated according the following expression, and 
20 output. 

[Expression 53] 
x = t • d 

M J ' MM **M J 

The likelihood calculating unit 60 
calculates errors using x M and x MJ (Expression 54) 
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according to Expression 55. 

[Expression 54] 

[Expression 55] 

After performing the error calculations, 
the likelihood calculating unit 60 receives the 
generated transmission symbol sequence corresponding 
to the transmitting antenna 11m/ receives 
corresponding error signals, holds S M transmission 
symbols having smaller numbers of errors 
corresponding to the transmitting antenna 11 M as 
described by Expression 56, and holds the 
corresponding errors as expressed by Expression 57. 
[Expression 56] 

[Expression 57] 

Then, the transmission symbol candidate 
generating unit 58 generates a transmission symbol 
as described by Expression 58 for the transmitting 
antenna ll M _i, and outputs. 

[Expression 58] 

d M-U 

The conversion signal replica generating 
unit 59 receives 

the R matrix, 

the S M transmission symbol sequences 
corresponding to the transmitting antenna 11m, and 

the transmission symbol corresponding to 
the transmitting antenna ll M -i 

for generating conversion signal replicas 



as described by Expression 59 according to 
Expression 60, and outputs the conversion signal 
replicas . 

[Expression 59] 
/\ 

X M-14 M -U M 

[Expression 60] 

x = r • d -h r • d 

Subsequently, the likelihood calculating 
unit 60 performs error calculations according to the 
following expression . 

[Expression 61] 

The likelihood calculating unit 60 holds 
Sm-i combinations of the transmission symbol 
sequences having smaller numbers of errors 
corresponding to the transmitting antenna 11 M and 
the transmitting antenna ll M -i as described by 
Expression 62 . 

[Expression 62] 

The likelihood calculating unit 60 holds 
errors at that time as described by Expression 63. 
[Expression 63] 

Similarly, the transmission symbol 
candidate generating unit 58 generates and outputs 
the transmission symbols d m ,i corresponding to the 
transmitting antenna ll m . The conversion signal 
replica generating unit 59 receives 

the R matrix, 

the transmission symbol corresponding to 
the transmitting antenna llm, and 
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Sm+i transmission symbol sequences 
corresponding to the transmitting antennas from ll m +i 
to 11„ 

for generating conversion signal replicas 
as described by Expression 64 according to 
Expression 65; and outputs the conversion signal 
replicas . 

[Expression 64] 
X 

[Expression 65] 

y — P * d ~\~ P * d ~h • ~f" v 1 

mA m mm mJ m m(m+l) m.i m+ i 'mi 

The likelihood calculating unit 60 
performs error calculations according to the 
following expression . 

[Expression 66] 



m.i m ,....ia/-1 J\f 



x — X 



m >'m ''m+l -~"*M 



+ e 



The likelihood calculating unit 60 holds 
S m combinations of the transmission symbol sequences 
corresponding to the transmitting antennas from ll m 
20 to 11 M as described by Expression 67, the 

combinations having smaller numbers of errors. 



[Expression 67] 



d 



'm(S m ). ■+ 1 4fl(S JI ). - ^.(^(5, 



.25 Further, the likelihood calculating unit 

60 holds the errors at that time as described by 
Expression 68 . 

[Expression 68] 



mJ m 0U m +\0) iM 0> mJ m (S m U m+ l(S m ) ft, (*«> 

30 By repeating the operations described 

above, CxS 2 combinations as described by Expression 
69 of the transmission symbol sequences 
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corresponding to all the transmitting antennas are 
obtained . 

[Expression 69] 

Further by repeating the operations 
described above, the errors at that time as 
described by Expression 70 are obtained. 

[Expression 70] 

Here, j is an integer equal to 1 or 
greater up to CxS 2 . 

The transmission symbol sequence 
estimating units 61i through 61 m receive surviving 
transmission symbols of the corresponding 
transmitting antennas, receive corresponding error 
signals, determine minimum errors, and estimate 
transmission symbol sequences that give the 
determined minimum errors . 

When the transmission symbol sequence 
estimating units 61i through 61 m carry out hard 
decisions, they are performed as follows. 

[Expression 71 ] 

d m =d m , m(J) where h(J),i 2 (j i M (j) provide- 
minimum of e mA(JXij{J) _ Ju{/) 

When the transmission symbol sequence 
estimating units 61i through 61 m carry out soft 
decisions, they are performed as follows. 

First, a bit sequence of soft decision i 
obtained according to the following expression. 

[Expression 72] 



u mj = tanh 
Here , 



mJ 
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[Expression 73] 
A., 

is a logarithmic likelihood ratio of a bit 
"i" of the transmitting antenna ll m , and is 
expressed as follows. 

[Expression 74] 

i 



Here, e m i n ,v is the smallest of 
[Expression 75] 

M0')^(/) (A/0) 

where the i-th bit is "V" , and a 2 (n) is 
noise power . 

The soft decision symbol 
[Expression 76] 



is estimated as shown by Expression 78, 
using Expression 77. 

[Expression 77] 

[Expression 78] 
Providing: d m =x m +j -y m 

In the case of QPSK J ! m " 

[y m = "m.2 

In the case of 16QAM < 



= iS" m4(2 ~" m ' 2) 



l) 

Finally, the multipath receiving signal 



demodulators 62i through 62 m receive the transmission 
symbol sequences that are estimated as described 
above, the transmission symbol sequences being 
described by Expression 79. 

[Expression 79] 



Then, the multipath receiving signal 
demodulators 62i through 62 ra multiply the received 
transmission symbol sequence as described by 
Expression 80, the spreading signal, and the channel 
coefficients as shown by Expression 82 such that 
received signal sequences, as described by 
Expression 81, for every reception path between the 
transmitting antennas and the receiving antennas are 
estimated and output. 

[Expression 80] 



[Expression 81] 

[Expression 82] 
L,(0=A>)-i(/> cit-r,) 

According to the multipath receiving 
signal demodulator 50, the number of times of error 

calculations is greatly reduced to C + ^ CS m ^ as 

described by Expression 83 in comparison with C M 
times required of the multipath receiving signal 
demodulator 40 that is shown by Fig. 3, although 
some degradation in demodulating the received signal 
sequence has to be tolerated. 

[Expression 83] 



Accordingly, complexity of the receiving 
apparatus is greatly reduced. 

Fig. 5 is a block diagram showing a 
multipath receiving signal demodulator 70 that is a 
fourth embodiment of the multipath receiving signal 
demodulator. According to the fourth embodiment, the 
multipath receiving signal demodulator uses the MLD 
algorithm as the primary demodulation method, 
wherein each path is individually processed. 

The multipath receiving signal demodulator 
70 includes a channel coefficient estimating unit 71, 
NxL reverse spreading units 72 through 75, L sets of 
QR factorizing units 76 and 77, L sets of Q H 
operating units 78 and 79 , a transmission symbol 
candidate generating unit 80, a conversion signal 
replica generating unit 81, a likelihood calculating 
unit 82, M transmission symbol sequence estimating 
units 83i through 83 m , and M multipath receiving 
signal demodulators 84i through 84 m . 

According to the multipath receiving 
signal demodulator 70 of the fourth embodiment, the 
channel coefficient estimating unit 71 receives a 
signal r n (t) received by the N receiving antennas 
21i through 21 n , and estimates the channel 
coefficients h m , n ,< of each path between the receiving 
antennas 21i through 21 n and the transmitting 
antennas Hi through ll m . 

Next, the reverse spreading units 72 
through 75, which number NxL, perform reverse 
spreading of the received signal r n (t) with the same 
spreading signal as used at the time of transmission 
at the receiving timing of each path, and NxL 
reverse-spread signals z n ,* one for each path of the 
receiving antennas . are obtained. 

Next, an f-th QR factorizing unit 
generates channel matrices that include channel 
coefficients of the *-th path, the number of the 
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channel matrices being equal to the number of the 
paths, and QR factorization of the channel matrices 
is performed such that a Q matrix and a R matrix are 
output . 

[Expression 84] 

\u Jh,\J ^3,1J ^4,U 

Channel Matrix H,= ^ ^ ^ ^ 

"l,3,/ "2.3J *h,3J ^4,3J 
Jh,4J klAJ ^3.4,/ h*AJ 



(where M = 4, N = 4) 



A Q< matrix output from the QR factorizing 
- unit is a unitary matrix having N rows x M columns, 
10 and suffices for Q* H Q<=I . Further, a R* matrix is an 
upper triangular matrix having M rows x M columns . 

Accordingly, an operation carried out by 
the *-th Q H operating unit can be described by 
Expression 85. 
15 [Expression 85] 



= Q, H z,= 



<ilU 



<ll2J <Jl.3J 















q*,*j_ 





Q, H (H,£> + N) = Q, H (Q,R,D + N) = R,jD + Q*N = 



0 
0 
0 



nzj 


nzj 




d x 


r 72J 


r 23J 


r 24J 


d 2 


0 






d> 


0 


0 




A. 



where M = 4, N = 4 



The transmission symbol candidate 
generating unit 80 generates the transmission 
20 symbols d M ,i corresponding to the transmitting 

antenna 11 M , and outputs them to the conversion 
signal replica generating unit 81. The conversion 
signal replica generating unit 81 receives the 
transmission symbols d M ,i, receives the R< matrix, 



- 38 - 



and generates conversion signal replicas as 
described by Expression 86 according to Expression 
87, and output the conversion signal replicas. 
[Expression 86] 

[Expression 87] 

The likelihood calculating unit 82 
calculates errors according to Expression 89 using 

x M ,< and x MJJ (Expression 88) . 

[Expression 88] 

[Expression 89] 

/«=1 

After performing the error calculation, 
the likelihood calculating unit 82 receives the 
generated transmission symbol sequences 
corresponding to the transmitting antenna 11 M / and 
the corresponding error signals. Then, the 
likelihood calculating unit 82 holds S M transmission 
symbols as described by Expression 90 having smaller 
numbers of errors corresponding to the transmitting 
antenna 11 M . 

[Expression 90] 

Further, the likelihood calculating unit 
82 holds errors at that time as described by 
Expression 9 1 . 

[Expression 91] 
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Next f the transmission symbol candidate 
generating unit 80 generates a transmission symbol 
as described by Expression 92 corresponding to the 
transmitting antenna ll M -i and outputs the 
transmission symbol . 

[Expression 92] 

The conversion signal replica generating 
unit 81 receives 

the R matrix, 

the S M transmission symbol sequences 
corresponding to the transmitting antenna 11 M , and 

the transmission symbol corresponding to 
the transmitting antenna ll M -i for generating 
conversion signal replicas as described by 
Expression 93, according to Expression 94, and 
outputs the conversion signal replicas. 

[Expression 93] 

[Expression 94] 

x =r • d +P • d 

Subsequently, the likelihood calculating 
unit 82 performs error calculations by the following 
expression . 

[Expression 95] 

e hf -\.i M _1 ,i M ~" ^ \ X M -1 ~~ *M -\J M _i .i^ | **" e M J M 

Then, the likelihood calculating unit 82 
holds Sm-i combinations as described by Expression 96 
of the transmission symbol sequences corresponding 
to the transmitting antenna 11 M and the transmitting 
antenna ll M -i, the combinations having smaller 
t numbers of errors. 

[Expression 96] 
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Further, the likelihood calculating unit 
82 holds the errors at that time as described by 
Expression 9 7 . 
5 [Expression 97] 

^-U w .,(l)>/ M (l)^^-U w . 1 (5 M . l )w M (5 w . 1 ) 

Similarly, the transmission symbol 
candidate generating unit 80 generates and outputs 
the transmission symbols d m ,i to the transmitting 
10 antenna ll m - The conversion signal replica 
generating unit 81 receives 

the R matrix, 

the transmission symbol corresponding to 
the transmitting antenna ll m , and 
15 S m +i transmission symbol sequences 

corresponding to the transmitting antennas from ll m +i 
to 11 M . 

Then, the conversion signal replica 

generating unit 81 generates conversion signal 

20 replicas as described by Expression 98 according to 

Expression 99 . 

[Expression 98] 

/\ 
X 

[Expression 99] 

The likelihood calculating unit 82 
performs error calculations according to the 
following expression . 

[Expression 100] 



30 



2 



The likelihood calculating unit 82 holds 
combinations as described by Formula 101 of the 



transmission symbol sequences corresponding to the 
transmitting antennas from ll m to 11 M , the 
combinations having smaller numbers of errors. 
[Expression 101] 

Further, the likelihood calculating unit 
82 holds the errors at that time as described by 
Expression 102. 

[Expression 102] 

m >'m (l)^m+l OX-*/ (1) m,i m {S m ),/ m+i {S m UJ M (S m ) 

By repeating the calculations described 
above, CxS2 combinations as described by Expression 
103 of the transmission symbol sequences 
corresponding to all the transmitting antennas are 
obtained . 

[Expression 103] 



Further, by repeating the calculations 
described above, the errors at that time as 
described by Expressing 104 are obtained. 

[Expression 104] 

where j is an integer equal to or greater 
than 1 up to CxS2 . 

The transmission symbol sequence 
estimating units 83 x through 83 m receive surviving 
transmission symbols and the corresponding error 
signals corresponding to each transmitting antenna, 
determine minimum errors, and estimate transmission 
symbol sequences that give the minimum errors. 

When the transmission symbol sequence 
estimating units 83i through 83 m carry out hard 
decisions, they are made as follows. 
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[Expression 105] 

d m = d mJm w where ). * 2 (/ ).-. ^ 0' ) Provide 
minimum e mMjWj) 

When the transmission symbol sequence 
estimating units 83i through 83 m carry out soft 
decisions, they are made as follows. First, a bit 
sequence of soft decision is obtained according to 
the following expression. 

[Expression 106] 



f 7 \ 



u mj = tanh 



10 Here, 

[Expression 107] 



is a logarithmic likelihood ratio of the 
bit "i" of the transmitting antenna ll m , and is 
15 described as follows. 

[Expression 108] 

i 



where e min ,v is the minimum of Expression 



109, 

20 [Expression 109] 



m ,q 0 ).r2 0 ..Jm 0 > r 



and ct 2 (n) expresses noise power. 
The soft decision symbol as described by 
Expression 110 is estimated as described by 
25 Expression 112 using Expression 111. 

[ Expression 110] 
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[Expression 111] 

[Expression 112] 
Providing :d m =x m +j -y m 



In the case of QPSK< 



= U m,2 



In the case of 16QAM 



5 Finally, the multipath receiving signal 

demodulators 84i through 84 m receive the 
transmission symbol sequences as described by 
Expression 113, the transmission symbol sequences 
being estimated as described above. 
10 [Expression 113] 



The transmission symbol sequences as 
described by Expression 114 are multiplied by the 
spreading signal and the channel coefficients so 

15 that received signal sequences as described by 

Expression 115 are obtained for every reception path 
from each transmitting antenna to every receiving 
antenna are obtained as described by Expression 116, 
and the received signal sequences are output. 

20 [Expression 114] 

d 

m 

[Expression 115] 
[Expression 116] 

25 LAO=L^>d m (ty c (t-r,) 
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According to the multipath receiving 
signal demodulator 70, similar to the case with the 
multipath receiving signal modulator 50 shown by Fig. 
4, the number of times of error calculations is* 

5 greatly reduced to C + ^GS^ as described by 

m«l 

Expression 117 in comparison with C M times required 
of the multipath receiving signal demodulator 40 
that is shown by Fig. 3, although some degradation 
in demodulating the received signal sequence has to 
10 be tolerated. 

[Expression 117] 

Fig. 6 is a block diagram of the channel 
coefficient estimating unit according to an 

15 embodiment of the present invention. Figs. 7A and 7B 
show frame structure examples of the signal 
transmitted from the transmitting apparatus when the 
channel coefficient estimating unit is used. 

First, the signal transmitted from the 

20. transmitting apparatus is described referring to 

Figs. 7A and 7B. According to the present embodiment, 
pilot symbols (indicated by hashed lines) , each 
being 4 symbols long, and each being different, are 
periodically inserted into data symbols of the 

25 - signals transmitted from the transmitting antennas 
(e.g., transmitting antennas 1 and 2) . In the 
example shown by Fig. 7A, pilot symbol patterns for 
the transmitting antennas are orthogonally arranged. 

Further, in the example shown by Fig. 7B, 

30 spreading signals (CI, C2) used for spreading the 
pilot symbols for the transmitting antennas are 
orthogonally arranged . 

By orthogonally arranging the pilot 
signals between the transmitting antennas using one 

35 of the methods shown by Figs. 7A and 7B, highly 



accurate channel estimation is attained. The 
orthogonal symbol patterns and the orthogonal 
spreading signals can be generated by using a Walsh 
sequence of the same length, for example, as the 
number of the pilot symbols, and the same length as 
the spreading signals for the pilot symbols. In the 
following description, a pilot signal sequence of a 
transmitting antenna m is expressed by p m (n), where 
n is a chip number. 

With reference to Fig. 6, the 
configuration of a channel coefficient estimating 
unit 100 according to the present embodiment is 
described. The channel coefficient estimating unit 
100 is for estimating a channel coefficient of each 
path between the transmitting antenna m and a 
receiving antenna n. Here, an example is shown 
wherein the channel coefficient estimating unit 100 
is for 3 transmitting antennas and 4 receiving 
antennas. That is, the channel coefficient 
estimating unit 100 includes 3x4 correlation 
detecting units 101 through 104, and pilot signal 
replica generating units 111 through 114. In 
addition, in this example, each of the correlation 
detecting units and the pilot signal replica 
generating units includes two or more components , 
which may be identified by a suffixed serial number 
when necessary. 

First, operations of estimating a channel 
coefficient hi,i,i between the transmitting antenna 1 
and the receiving antenna 1 are described with 
reference to Fig. 6. 

A signal rl is received by the receiving 
antenna 1, and is provided to the correlation 
detecting unit 101. Further, the pilot signal 
replica generating unit 111 generates a pilot symbol 
sequence pi corresponding to the transmitting 
antenna 1, which pi is provided to the correlation 



- 46 - 



detecting unit 101. 

The correlation detecting unit 101 
estimates the channel coefficient hi.i.i between the 
transmitting antenna 1 and the receiving antenna 1 
5 according to Expression 118, wherein 

the received signal rl is multiplied by a 
complex-conjugate value of the pilot symbol sequence 
pi of the transmitting antenna 1 in consideration of 
the receiving timing of a path t , and 
10 the obtained product is averaged over the 

4 pilot symbol periods. 

[Expression 118] 

Here, r x (n) represents the received signal 

15 ri when a pilot symbol n is received. An 

implementation is possible wherein hi.i.i is 
estimated by weight-averaging channel coefficients 
obtained in two or more pilot symbol periods. 

Similarly, a second step correlation 

20 detecting unit 102 receives the signal rl, and a 

pilot symbol sequence p 3 of the transmitting antenna 
3 generated by a second step pilot symbol replica 
generating unit 112, and a channel coefficient hi, 3 ,< 
is estimated and output. 

25 Further, a received signal r4 and the 

pilot symbol sequence pi of the transmitting antenna 
1 generated by the pilot symbol replica generating 
unit 113 are provided .to the correlation detecting 
unit 103, and correlation is obtained such that a 

30 channel coefficient h 4 ,i,i is estimated and output. 

By repeating the above operation, the 
channel coefficient of all the paths between the 3 
transmitting antenna and the 4 receiving antennas 
are estimated. According to the example described 

35 above, the pilot symbol is multiplexed by the data 

symbol in the time domain; nevertheless, the channel 
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coefficient estimation can be obtained by the same 
method in the case where code multiplexing is 
employed . 

Fig. 8 is a block diagram of a receiving 
5 apparatus 200 that is a second embodiment of the 
receiving apparatus of the present invention. As 
shown in Fig. 8, the receiving apparatus 200 
includes multipath receiving signal demodulators 211 
through 213 that are serially connected through 

10 multipath interference canceling units 221i, 221 2/ 
222i, 222 2 , 223i, and 223 2 (i.e., this example is a 
3-stage configuration) , and a demodulator 231 that 
is arranged at the last stage. To the multipath 
receiving signal demodulator 211 , serving as the 

15 first stage, any one of the multipath receiving 

signal demodulators described with reference to Figs. 
2 through 5 can be applied. 

Further, any one of the multipath 
receiving signal demodulators described with 

20 reference to Figs. 2 through 5 can be applied to the 
multipath receiving signal demodulators 212, 213, 
and so on, as applicable. Here, signals that are 
input to the p-th multipath receiving signal 
demodulator (p is an integer equal to 2 or greater) 

25 are described by Expression 119. 

[Expression 119] 

This can be described as Expression 121, 
using the received signals, and output signals 
30 (Expression 120) of the (p-l)-th multipath receiving 
signal demodulator . 

[Expression 120] 

Specifically, the signals that are input 
35 to the p-th multipath receiving signal demodulator 



are generated by the preceding multipath 
interference canceling units according to the 
following expression . 

[Expression 121] 

NxL received signals after multipath 
interference cancellation obtained according to 
Expression 121 are described by Expression 122. 

[Expression 122] 

The NxL received signals are provided to a 
reverse spreading unit (not illustrated) 
corresponding to the receiving antenna 21 n and the 
path i . 

Further, according to the configuration 
shown by Fig. 8, a channel coefficient estimating 
unit (not illustrated) of the p-th stage multipath 
receiving signal demodulator may use a received 
signal after multipath interference cancellation as 
described by Expression 123 instead of the received 
signal r n (t) when estimating the coefficient hi.i.i, 
which results in a more accurate channel coefficient 
estimation . 

[Expression 123] 

(It is possible to use the channel 
coefficient of the (p-l)-th multipath receiving 
signal demodulator for simplification of the 
configuration . ) 

As described above, by preparing two or 
more stages of the multipath receiving signal 
demodulators according to the present embodiment, 
the latter stage multipath receiving signal 



demodulators can use the received signal after 
multipath interference cancellation, which enables 
highly accurate channel estimation and transmission 
symbol sequence estimation. As a result, the 
received signal sequence as described by Expression 
124 is accurately estimated for every reception path 
between the transmitting antennas and the receiving 
antennas . 

[Expression 124] 

Fig. 9 is a block diagram of a demodulator 
300 that is a first embodiment of the demodulator 
applicable to the receiving apparatus according to 
the embodiments of the present invention. The 
present embodiment uses MLD for the demodulator as 
the demodulating algorithm. 

The demodulator 300 includes a channel 
coefficient estimating unit 311 , NxL reverse 
spreading units 312 through 315, a transmission 
symbol candidate generating unit 316, received 
signal replica generating unit 317, a likelihood 
calculating unit 318, and a transmitting sequence 
estimating unit 319. Where a certain unit includes 
two or more components, they are identified with a 
suffixed serial number. 

An operation of the demodulator 300 is 

described . 

The demodulator 300 is provided with NxL 
received signals after multipath interference 
cancellation as described by Expression 125. 

[Expression 12 5] 

The channel coefficient estimating unit 
311 estimates the channel coefficients h m , n ,< of each 
path between the receiving antennas and the 
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transmitting antennas (here, it is also possible to 
use a channel coefficient of a preceding multipath 
receiving signal demodulator for simplification of 
the configuration) . 

Further, NxL received signals after 
multipath interference cancellation as described by- 
Expression 126 are provided to a reverse spreading 
unit corresponding to a receiving antenna 21 n and 
the path I such that reverse-spread signals z'n.a 
are obtained. 

[Expression 126] 



The transmission symbol candidate 
generating unit 316 generates and outputs the 
transmission symbols d m ,i corresponding to each 
antenna. The received signal replica generating unit 
317 receives the transmission symbols and the 
channel coefficients, and generates received signal 
replicas as described by Expression 127 according to 
Expression 128, and outputs the received signal 
replicas to the likelihood calculating unit 318. 

[Expression 127] 

Z 

nJA ,i2 *M 

[Expression 128] 

The likelihood calculating unit 318 
receives, the reverse-spread signals z' n ,i from the 
reverse spreading units 312 through 315, and the 
received signal replicas as described by Expression 
129 . 

[Expression 129] 
Z 

Then, the likelihood calculating unit 318 



performs error calculations according to the 
following expression . 

[Expression 130] 



e =YYk,-£ 



The generated transmission symbols 
corresponding to each transmitting antenna, and the 
corresponding error signals are provided to the 
transmitting sequence estimating unit 319 that 
outputs a likelihood A-i of a bit transmitted by the 
transmission symbol sequence. Here, the bit 
likelihood based on the error signal may be 
calculated by any available method. 

The bit likelihood described above is 
provided to a channel decoder (for example, a turbo 
decoder), etc., so that an information bit sequence 
is finally obtained. 

Fig. 10 is a block diagram of a 
demodulator 400 that is a second embodiment of the 
demodulator of the present invention. The second 
embodiment uses MLD as the demodulating algorithm, 
wherein paths are grouped into a block for 
processing. 

The demodulator 400 includes a channel 
coefficient estimating unit 411, NxL reverse 
spreading units 412 through 415, a QR factorizing 
unit 416, a Q H operating unit 417, a transmission 
symbol candidate generating unit 418, a conversion 
signal replica generating unit 419, a likelihood 
calculating unit 420, and a transmitting sequence 
estimating unit 421. Where a certain unit includes 
two or more components, they are identified with a 
suffixed serial number. 

Operations of the demodulator 400 are 
described below. 

The demodulator 400 receives NxL signals 
as described by Expression 131 after multipath 
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10 



15 



20 



interference cancellation. 

[Expression 131] 



The channel coefficients h m , n ,< of each 
path of between the receiving antennas and the 
transmitting antennas are estimated by the channel 
coefficient estimating unit 411. (Here, it is also 
possible similarly to use a channel coefficient in a 
preceding multipath receiving signal demodulator for 
simplification of the configuration) . 

Furthermore, the NxL received signals 
after multipath interference cancellation 

[Expression 132] 



do 



are provided to the reverse spreading unit 
corresponding to the receiving antenna 21 n and path 
t, and reverse-spread signals z' n , fl are obtained. 

Next, the QR factorizing unit 416 
generates the following channel matrix (Expression 
133) containing channel coefficients, performs QR 
factorization of the channel matrix, and outputs a Q 



matrix and a R matrix to the Q 
[Expression 133] 



operating unit 417. 



Channel Matrix H = 





^2,1,1 


^3.1.1 


^4,1,1 


^1,2,1 


^.2.1 


^3.2,1 


^4,2,1 


^.3,1 


^,3,1 


^3,3.1 


^4,3.1 




^2,4,1 


^3.4.1 


^4,4,1 


K\a 


^2,1,2 


^3,1.2 


^4,1,2 


^1.2.2 


^.2,2 


^3,2,2 


^4.2,2 


^1.3.2 


^.3,2 


KX2 


^4,3,2 


^1.4,2 


^2,4.2 


^3,4,2 


^4,4,2 


4,N 


= 4,L = 


2) 





25 



H = QR 

The Q matrix output from the QR 
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factorizing unit 417 is a unitary matrix having 
(NxL) rows x M columns, and suffices Q H Q=I . Further, 
the R matrix is an upper triangular matrix having M 
rows x M columns. Accordingly, an operation carried 
out by the Q H operating unit 417 can be described as 
follows. 

[Expression 134] 



#1,1,1 #1,2,1 #1,3,1 #1.4,1 

* • * * 

#2.1,1 #2,2,1 #2,3,1 #2,4,1 

* • * * 
#3,1,1 #3,2,1 #3.3,1 #3,4,1 

#1.1.1 #4.2,1 #1,3,1 #4,4,1 









X = 


*2 


= Q H Z = 










_*4_ 





Q H (HZ) + N) = Q H (QRD + N) = R£> + Q H N = 













z \.\ 












z 2 .i 


#U.2 


#1*2,2 


#1*3.2 


#1.4.2 


z \\ 


#2.1.2 


* 

#2.2,2 


#2.3.2 


#2.4.2 


z *.i 


#3.1.2 


* 

#3.2.2 


#3.3.2 


#3.4.2 


Z 1.2 


#1.1.2 


#4.2,2 


#4.3.2 


• 

#4.4,2 _ 


Z 2,2 












Z 3.2 












Z 4.2 
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"^11 


r u r i3 


^14 " 
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r 22 r 23 


>24 




+ 




0 


0 r 33 


>34 


d 3 
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0 0 


>44_ 


A. 




< 



(In this example, M = 4, N = 4, L = 2) 

The transmission symbol candidate 
10 generating unit 418 generates and outputs the 
transmission symbols d M ,j corresponding to. the 
transmitting antenna 11 M . The conversion signal 
replica generating unit 419 receives the 
transmission symbols d M ,j and the R matrix. Then, the 
15 conversion signal replica generating unit 419 

generates conversion signal replicas as described by 
Expression 135 according to Expression 136, and 
outputs the conversion signal replicas to the 
likelihood calculating unit 420. 
20 [Expression 135] 



"•J 



[Expression 136] 
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r = r • u 

M J AO/ A/ J 

The likelihood calculating unit 420 

calculates errors between x' M and x MJ (Expression 

137) according to Expression 138. 
[Expression 137] 



[Expression 138] 
I r _ - I 2 

^MJ \ X M X MJ\ 



The generated transmission symbol 
10 sequences corresponding to the transmitting antenna 
11 M and corresponding error signals are provided to 
the likelihood calculating unit 420 for generating 
S M transmission symbols as described by Expression 
139 having smaller numbers of errors corresponding 
15 to the transmitting antenna 11 M . 

[Expression 139] 

The likelihood calculating unit 420 
further calculates errors at that time as shown by 
20 Expression 140. 

[Expression 140] 

The transmission symbol candidate 
generating unit 418 generates the transmission 
25 symbols as described by Expression 141 corresponding 
to the transmitting antenna ll M -i, and outputs the 
transmission symbols . 

[Expression 141] 

3 0 The conversion signal replica generating 
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unit 419 receives 

the R matrix, 

the S M transmission symbol sequences 
corresponding to the transmitting antenna 11 M and, 
5 the transmission symbols corresponding to 

the transmitting antenna ll M -i for generating 
conversion signal replicas as described by 
Expression 142 using Expression 143, and the 
conversion signal replicas are output. 
10 [Expression 142] 

[Expression 143] 

x — p • H +p • d 

Subsequently, the likelihood calculating 
15 . unit 420 performs error calculations according to 
the following expression. 

[Expression 144] 



e M-iJ M ^J Ki \ X M-l X M-Um-iJm I ~^ e MJM 

The likelihood calculating unit 420 holds 
20 S M -i combinations of the transmission symbol 

sequences corresponding the transmitting antenna 11 M 
and the transmitting antenna ll M -i, the combinations 
having smaller numbers of errors. 
[Expression 145] 

Further, the likelihood calculating unit 
420 holds the errors at that time as described by 
Expression 146. 

[Expression 146] 

3 0 ^w-i^-i(n^o) e u - 1 -i -i >-w -i > 

Similarly, the transmission symbol 
candidate generating unit 418 generates and outputs 



the transmission symbols d m ,i corresponding to the 
transmitting antenna ll m . The conversion signal 
replica generating unit 419 receives 

the S m +i transmission symbol sequences 
corresponding to the transmitting antennas from ll m +i 
to 11 M , 

the transmission symbol corresponding to 
the transmitting antenna ll m , and 

the R matrix for generating conversion 
signal replicas as described by Expression 147 
according to Expression 148 , and output the 
conversion signal replicas. 

[Expression 147] 

[Expression 148] 

m,i m .i m +i...J M mm m.i m m(m+\) m^Wl ™M m **M 

The likelihood calculating unit 420 
performs error calculations according to the 
following expression. 

[Expression 149] 

f ~ I 2 , 

£ = X — X + B 

rn.i m i M _! j M m m ,i m i M | m+U m +l -'m + 2 

The likelihood calculating unit 420 holds 
S m combinations as described by Expression 150 of 
the transmission symbol sequences corresponding to 
the transmitting antennas from ll m to 11 M , the 
combinations having smaller numbers of errors. 

[Expression 150] 

Further, the likelihood calculating unit 
420 holds the errors at that time as described by 
Expression 151. 

[Expression 151] 
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By repeating the operations described 
above, CxS 2 combinations of the transmission symbol 
sequences corresponding to all the transmitting 
5 antennas as described by Expression 152 are obtained, 
[Expression 152] 

{d^ .o ) ' .o > 9 " ' • ' d MJM , 0 ) } 

Further, the errors as described by 
Expression 153 at that time are obtained. 
10 [Expression 153] 

m,f| 

The transmitting sequence estimating unit 
421 receives the surviving transmission symbols 
corresponding to each transmitting antenna and the 

15 corresponding error signals, and outputs a 
likelihood of a bit transmitted by the 
transmission symbol sequence. Here, the bit 
likelihood based on the error signal can be 
calculated using any available method. 

20 The bit likelihood described is provided 

to a channel decoder (for example, a turbo decoder) 
etc., and an information bit sequence is finally 
obtained . 

According to the demodulator 400, the 
25 number of times of error calculations is greatly 

M -I 

reduced to C + CS^ as described by Expression 154 

in comparison with C M times required of the 
demodulator 300 that is shown by Fig. 9, although 
some degradation in the bit likelihood estimation 
30 has to be tolerated. 

[Expression 154] 
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In this way, the complexity of the 
receiving apparatus can be greatly decreased. 

Fig. 11 shows a demodulator 500 that is a 
third embodiment of the demodulator of the present 
invention. According to the demodulator 500 of the 
third embodiment, each path is processed by a 
demodulating algorithm that employs MLD . 

The demodulator 500 includes a channel 
coefficient estimating unit 511, NxL reverse 
spreading units 512 through 515, L QR factorizing 
units 516 and 517, L Q H operating units 518 and 519, 
a transmission symbol candidate generating unit 520, 
a conversion signal replica generating units 521, a 
likelihood calculating units 522, and a transmitting 
sequence estimating unit 523. Where a certain unit 
includes two or more components, they are identified 
with a suffixed serial number. 

Operations of the demodulator 500 are 
described below. 

The demodulator 500 receives NxL input 
signals that are the received signals after 
multipath interference cancellation as described by 
Expressed 155 . 

[Expression 155] 

The channel coefficient estimating unit 
511 estimates the channel coefficients h m , n , fi of each 
path between the receiving antennas and the 
transmitting antennas. (Here, it is possible to use 
a channel coefficient of a preceding multipath 
receiving signal demodulator for simplification of 
the configuration) . 

Furthermore, the NxL received signals 
after multipath interference cancellation as 
described by Expression 156 are provided to the 
reverse spreading unit corresponding to the 



receiving antenna 21 n and the path I so that the 
reverse-spread signals z' n ,fl are obtained. 
[Expression 156] 

Next, the f-th QR factorizing unit 
generates the following channel matrix containing 
channel coefficients of the *-th path for each path, 
performs QR factorization of the channel matrix so 
that the Q matrix and the R matrix are output. 

[Expression 157] 

Channel Matrix H t 

(Example where M 
H, = Q,R, 

A Qt matrix output by the QR factorizing 
unit is a unitary matrix having N rows x M columns, 
and suffices Q* H Q*=I . Here, H expresses a conjugate 
complex transposition, and "I" expresses a unit 
matrix. Further, a R< matrix turns into an upper 
triangular matrix having M rows x M columns. 

Accordingly, the operation carried out by 
the *-th Q H operating unit can be described as 
follows . 

[Expression 158] 
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(Example where M = 4, N = 4) 

The transmission symbol candidate 
generating unit 520 generates the transmission 
symbols d M ,i corresponding to the transmitting 
antenna 11 M , and outputs the transmission symbols 
d M ,i to the conversion signal replica generating unit 
521. The conversion signal replica generating unit 
521 receives the transmission symbols d M ,i and the R* 
matrix so that the conversion signal replicas as 
described by Expression 159 are generated according 
to Expression 160, and output. 

[Expression 159] 
/\ 

X MJ 

[Expression 160] 



X MJL4 r * 



20 



MJJ 

The likelihood calculating unit 522 
performs error calculations using x'm,* and the 
conversion signal replicas as described by 
Expression 161 according to Expression 162. 

[Expression 161] 



'MJJ 



[Expression 162] 



e„. = 



= Vlv' _ 



M JJ I 



After the error calculations, the 
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likelihood calculating unit 522 receives generated 
transmission symbol sequences corresponding to the 
transmitting antenna 11 M , and the corresponding 
error signals. The likelihood calculating unit 522 
5 holds S M transmission symbol sequences having 
smaller numbers of errors corresponding to the 
transmitting antenna 11 M as described by Expression 
163. 

[Expression 163] 

Further, the likelihood calculating unit 
522 holds the errors at that time as described by 
Expression 164 . 

[Expression 164] 

15 ^Mi(l)^Wj%) 

Next, the transmission symbol candidate 
generating unit 520 generates and outputs 
transmission symbols d M _i ri corresponding to the 
transmitting antenna ll M -i- 
20 The conversion signal replica generating 

unit 521 receives 

the R matrix, 

the S M transmission symbol sequences 
corresponding to the transmitting antenna 11 M and 

25 the transmission symbol corresponding to 

the transmitting antenna ll M -i for generating 
conversion signal replicas as described by 
Expression 165 according to Expression 166, and 
outputting the conversion signal replicas. 

30 [Expression 165] 

[Expression 166] 

y ~ r * d + r • d 

■^M-\.IJ M -\,i M '(A/-DCA/-D./ **M-U M -i '(Af-DM ./ M .i M 



Subsequently the likelihood calculating 
unit 522 performs error calculations according to 
the following expression. 

[Expression 167] 

The likelihood calculating unit 522 holds 
S M -i combinations as described by Expression 168 of 
the transmission symbol sequences corresponding the 
transmitting antenna 11 M and the transmitting 
antenna ll M -i having smaller numbers of errors. 

[Expression 168] 

Further, the likelihood calculating unit 
522 holds the errors at that time as described by 
Expression 169. 

[Expression 169] 

p P 

Similarly, the transmission symbol 
candidate generating unit 520 generates and outputs 
the transmission symbols d m ,i corresponding to the 
transmitting antenna ll m . The conversion signal 
replica generating unit 521 receives 

the R matrix, 

S m+ i transmission symbol sequences 
corresponding to the transmitting antennas m+1 to M, 
and 

the transmission symbol corresponding to 
the transmitting antenna ll ra for generating 
conversion signal replicas as described by 
Expression 170 according to Expression 171, and 
output . 

[Expression 170] 
X 



[Expression 171] 

v — r 9 d 4- y • rj -L -4- P * d 

The likelihood calculating unit 522 
performs error calculations according to the 
following expression . 

[Expression 172] 

e = V y' - x + e 

/el 

The likelihood calculating unit 522 holds 
. S m combinations of the transmission symbol sequences 
corresponding to the transmitting antennas ll m to 
11 M having smaller numbers of errors as described by 
Expression 173. 

[Expression 173] 

Further, the likelihood calculating unit 
5-22 holds the errors at that time as described by 
Expression 174. 

[Expression 174] 

& at 

mJ m (DJ m +\ 0) ft/ (1) m ,i m (S m )./ OT+ i (S m ) ft/ (S m ) 

By repeating the above operations, CxSz 
combinations of the transmission symbol sequences 
corresponding to all the transmitting antennas, and 
the errors at that time are obtained as described by 
Expressions 175 and 176, respectively. 

[Expression 175] 

[Expression 176] 
& 

m.q(y).i2 0) WO) 

The transmitting sequence estimating unit 
523 receives the surviving transmission symbols 
corresponding to each transmitting antenna obtained 



by the L likelihood calculating units, and the 
corresponding error signals; and outputs a 
likelihood transmitted by the transmission symbol 
sequence based on the sum of L error signals. Here, 
the bit likelihood based on the error signals can be 
calculated using any available method. 

The bit likelihood described above is 
provided to a channel decoder (for example, turbo 
decoder) , etc. , and, the information bit sequence is 
finally obtained. 

According to the demodulator 500 (Fig. 11) 
the number of times of error calculations is greatly 

reduced to C + ^ CS m+l as described by Expression 177 

mml 

in comparison with C M times required of the 
demodulator 300 that is shown by Fig. 9, although 
some degradation in the bit likelihood estimation 
has to be tolerated. 

[Expression 177] 

In this way, the complexity of the 
receiving apparatus is greatly reduced. 

Fig. 12 is a block diagram of a receiving 
apparatus 600 that is a third embodiment of the 
receiving apparatus according to. the present 
invention. The receiving apparatus 600 constitutes a 
radio communications system wherein a transmitting 
apparatus (not illustrated) performs serial parallel 
conversion of a transmitting data bit sequence into 
M transmitting sequences, which sequences are 
modulated with the same frequency and the same 
spreading code group, and are simultaneously 
transmitted by N CO d e multicode channels from M 
transmitting antennas . 

The receiving apparatus 600 includes 
multipath receiving signal demodulators 611 through 
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614 corresponding to each code channel (each of code 
channels 1 and 2 in this example) for a 
predetermined number of stages, demodulators 621 and 
622 corresponding to each code channel, and 

5 multipath interference canceling units 615 through 
618. The multipath interference canceling units 615 
through 618 are inserted between one and another of 
the multipath receiving signal demodulators 611 
through 614 , or between one of the multipath 

10 receiving signal demodulators 611 through 614 and 

one of the demodulators 621 and 622, as applicable. 

According to the present embodiment, the 
multipath receiving signal demodulators 611 and 612 
that represent the first stage of each code channel 

15 receive signals r n (t) that are received by 

corresponding N receiving antennas (N=2 in this 
example) , and output demodulated signal sequences as 
described by Expression 178 for each reception path 
from the transmitting antennas corresponding to the 

20 code channels. 

[Expression 178] 

i'::,st)=k^;dr(.t> c t (. t -T,) 

where k represents a code number that is 
limited by Expression 179. 
25 [Expression 179] 

The multipath interference canceling units 

615 and 616 (note that N=2 in this example) receive 
the signals r n (t), and demodulated signal sequences 

30 of all the code channels as described by Expression 
180 through corresponding adders 631 through 634. 
[Expression 180] 

The multipath interference canceling units 



615 and 616 then generate received signals after 
multipath interference cancellation as described by 
Expression 181 according to Expression 182. 
[Expression 181] 

[Expression 182] 

The multipath receiving signal 
demodulators 613 and 614 that are the second stage 
of the corresponding code channels receive the 
received signals after multipath interference 
cancellation as described by Expression 183. 

[Expression 183] 

<«. 

The operations described are repeated for 
subsequent stages as applicable, and final received 
signals after multipath interference cancellation as 
described by Expression 184 are generated. 

[Expression 184] 




The final received signals after multipath 
interference cancellation generated as above are 
provided to the corresponding demodulators 621 and 
622 of each code channel, and the demodulators 621 
and 622 output the likelihood of a bit being 
transmitted by the transmission symbol sequence in 
each code channel . 

Next, computer simulations of the 
embodiments of the present invention are described 
with reference to Fig. 13. Fig. 13 shows computer 
simulation results of throughput property (vertical 
axis) over average receiving Eb/N D per receiving 



antenna (signal energy to the noise power density 
per 1 bit of information) (horizontal axis) of 
various demodulating methods used by the receiving 
apparatus; namely, the conventional MLD , the 
conventional 2-dimensional MMSE , and the method of 
the present invention. 

The receiving apparatus applied to the 
computer simulation includes two stages of the 
multipath receiving signal demodulators, where the 
configuration of Fig. 2 is used for the first stage, 
and the configuration of Fig. 4 is used for the 
second. Further, as for the multipath, a two-path 
model is used where both paths provide the same 
average received power. 

The spreading rate is 16, and 15 code 
channels are code-multiplexed (effective spreading 
rate is 15/16) . In Fig. 13, a coding rate of channel 
coding using turbo coding is indicated by "R" , solid 
lines represent properties . according to the present 
invention (QR-MLD with two-stage MPIC) , dotted lines 
represent properties of the conventional MMSE, and X 
marks represent properties of the conventional MLD. 
Further, • and O marks represent the average 
receiving Eb/No per receiving antenna in the case of 
MIMO multiplexing (4 transmitting antenna , 4 
receiving antennas) using QPSK modulation; ■ and □ 
marks represent the average receiving Eb/No per 
receiving antenna in the case of MIMO multiplexing 
using 16QAM modulation. 

As shown by Fig. 13, average received 
power Eb/No required for obtaining a certain level 
of the throughput is greatly reduced by using the 
configuration of the receiving apparatus according 
to the embodiments of the present invention as 
compared with the case where the conventional signal 
separating method is used. That is, the receiving 
apparatus of the present invention requires less 



transmission power to obtain the same throughput as 
conventional apparatuses. In other words, the same 
transmitted power provides a greater throughput with 
the receiving apparatus according to the embodiments 
of the present invention. 

As described above, according to the 
present invention, when the MIMO multiplexing is 
used for simultaneous transmission of different data 
from two or more transmitting antennas in order to 
increase an information bit rate in a CDMA system, 
highly accurate signal separation is obtained, 
reducing the multipath interference. Consequently, 
the receiving bit error rate and receiving packet 
error rate in a multipath fading environment are 
greatly reduced, and the throughput (information bit 
transmission rate that provides no errors) is 
greatly raised. 

Although the multipath receiving signal 
demodulator has been described as using algorithms 
(for example, a program that describes the 
algorithm) such as MMSE and MLD , other algorithms 
such as QR-MLD algorithm can be more desirable for 
obtaining accurate signal separation by a realistic 
amount of calculations. 

Further, although the receiving apparatus 
has been described as being used in a radio 
communications system, the receiving apparatus can 
be applied to a base station of a mobile 
communications system. Further, the algorithm. of the 
present invention can be applied to a so-called 
software radio base station, to which the algorithm 
is selectively downloaded to a radio apparatus 
through a radio circuit or a wire circuit for 
changing properties of the radio apparatus . 

[Availability to industry] 

The present invention is applicable to a 
radio communications system that demodulates signals 
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using two or more antennas. 

The present application is based on 
Japanese Priority Application No. 2004-063197 filed 
on March 5, 2004 with the Japanese Patent Office, 
5 the entire contents of which are hereby incorporated 
by reference . 



